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Abstract. Dynamics of external cavity semiconductor lasers is known to be a 
complex and uncontrollable phenomenon. Due to the lack of experimental studies on 
the nature of the external cavity semiconductor lasers, there is a need to theoretically 
clarify laser dynamics. The stability of laser dynamics in the present paper, is analyzed 
through plotting the Lyapunov exponent spectra, bifurcation diagrams, phase portrait 
and electric field intensity time series. The analysis is preformed with respect to applied 
feedback phase C p , feedback strength r\ and the pump current of the laser. The main 
argument of the paper is to show that the laser dynamics can not be accounted for 
through simply a bifurcation diagram and single-control parameter. The comparison 
of the obtained results provides a very detailed picture of the qualitative changes in 
laser dynamics. 
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1. Introduction 

Semiconductor lasers are ubiquitous in the modern world. They are the principal source 
of coherent light in optical communications and ultra-fast optical processing, they are 
used in optical storage devices,in laser pointers, to cite only a few applications [11 EJ [7] . 
Semiconductor lasers provide one of the best physical systems for studying nonlinear 
dynamic phenomena. Many researchers have investigated the chaotic behavior of several 
laser systems [H El E] . 

Semiconductor laser has many advantages, such as small size, easy integration, 
compactness, low cost and convenience of operation. Therefore, they are preferable 
to any other types of lasers in the field of optical telecommunications. 
External cavity semiconductor lasers (ECSLs) are an integral part of high speed chaos 
based communication systems [El [9]. Hence, ECSLs have been a subject of extensive 
research [101 [TTJ [12] . Understanding the influence of delayed optical feedback on the 
behavior of ECSLs is of great relevance for technological applications. 
Different studies have been conducted to characterize or manipulate of ECSLs; long 
external cavity [T2J, [14] and short external cavity [TBI IT6l IT?] . Because of advantages in 
the short cavity regime [El [18], in this paper, we just focus on the short external cavity 
semiconductor lasers as a dynamical system. The main model in the understanding 
of the dynamics of these lasers is the well-known Lang-Kobayashi (LK) equations, 
consisting of two delay differential equations (DDEs) [21] for the complex electrical 
field E and the carrier number density N [T9J [20] . These equations closely describe the 
chaotic behavior of the physical system and confirm the presence of chaotic regimes [13] . 
A large part of these studies have been made possible using numerical continuation 
methods which allow one to find and follow numerical solution to the rate equations 
irrespective of their stability. 

Our main point here is that one is likely to miss important phenomena if one just 
considers a bifurcation diagram. This is why we show the behavior of the laser for a 
single value of the control parameter in different ways: by Lyapunov exponents and 
time series. This allows us to present a consistent overall picture of the dynamics at 
the same time. In this study parameters such as feedback phase, feedback strength, 
and pump current in the short external cavity are used as control parameters to obtain 
different dynamical regimes, including periodic and quasi periodic (QP) [22], regular 
pulse packages (RPP) [IB], and chaotic (CH) behavior [23] . 

2. The Laser model 

In the early 1980s, Lang-Kobayashi (LK) proposed model semiconductor lasers 
[241 125] . The LK equations are model equations that have been used extensively in the 
past to describe a semiconductor laser subject to feedback from an external cavity [25] 
where they used DDEs with the advantage that have an infinite-dimensional phase 
space [21] . For the (complex) electric field E and inversion N, we write the LK equations 
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as the dimensionless and compact set of equations 
dE 

— = (l + ia)NE + nE(t-T)e- lCp (1) 

T^j- = P- N -(1-2N) \E\ 2 . (2) 

The parameters in the above-mentioned equations describe the line width enhancement 
factor a, the feedback strength r] , the 27r-periodic feedback phase C p , the ratio between 
carrier, photon lifetime T and the pump current P. Equations (1), (2) describe a 
semiconductor laser with external optical feedback. In these equations, the time is 
normalized to the cavity photon lifetime ( lps) and T is the ratio of the carrier lifetime 
(Ins) to the photon lifetime [IB] . In this study the parameters are chosen not only to 
elucidate the dynamical structure but also correspond fairly well to the experimental 
conditions. The external round trip time r is also normalized to the photon lifetime. 
The remaining parameters, however, are held fixed at T = 1710, r = 70, a = 5.0 [15J. 
In this paper the stability of an electric field intensity | E p is studied versus Cp, rj , 
and P. Where all the parameters are easily accessible in experiments [T5"| 126]. 



3. Stability analysis 

3.1. Bifurcation diagrams 

Bifurcation means a qualitative change in the dynamical behavior of a system when a 
parameter of the system is varied. A bifurcation diagram provides a useful insight into 
the transition between different types of motion that can occur as one parameter of the 
system alters [27]. It enables one to study the behavior of the system on a wide range 
of an interested control parameter. In this paper the dynamical behavior of the system 
is studied through plotting the bifurcation diagrams of the | E | 2 versus C p , rj, and P 
as control parameters. This procedure continued by increasing the control parameters, 
and the new resulting points were plotted in the bifurcation diagram versus the new 
control parameter. 



3.2. Lyapunov exponent spectrum 

Lyapunov exponents and entropy measures, on the other hand can be considered as 
'dynamic' measures of attractor complexity and they are called 'time average' [27J. 
The Lyapunov exponent A is useful for distinguishing various orbits. Lyapunov 
Exponents quantify sensitivity of the system to initial conditions and give a measure of 
predictability. The Lyapunov exponent is a measure of the rate at which the trajectories 
separate one from another. A negative exponent implies that the orbits approach to 
a common fixed point. A zero exponent means that the orbits maintain their relative 
positions; they are on a stable attractor. Finally, a positive exponent implies that 
the orbits are on a chaotic attractor, so the presence of a positive Lyapunov exponent 
indicates chaos. The Lyapunov exponent is defined as follows: 
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Consider two nearest neighboring points in phase space at time and t, with distances 
of the points in the ith direction ||&Ej(0)||, and respectively. The Lyapunov 

exponent is then defined by the average growth rate A$ of the initial distance, 



Aj = lim - In J J _ WC'i (3) 
t^oo t || 5^(0)11 v ; 

The existence of a positive Lyapunov exponent is the indicator of chaos showing 
neighboring points with infinitesimal differences at the initial state abruptly separate 
from each other in the ith direction [28]. Using the algorithm of Wolf [29], the Lyapunov 
exponent was calculated versus a given control parameter. Then the value of the control 
parameter increased a little and the Lyapunov exponent was calculated for the new 
control parameter. By continuing this procedure Lyapunov exponent spectrum of the 
system was plotted versus the control parameter. 



4. Results and Discussions 



The absence of employing direct mathematical methods in study of ECSLs with respect 
to the variations of the parameters C p , rj, and P easily could be solved by consider- 
ing the bifurcation diagram and Lyapunov exponent spectrum. The output intensity 
(| E | 2 ) dynamics in the ECSLs, has been studied by considering both the variation of 
each parameter (rj, C P ,P) and the two other (77, C p , P) as a initial condition of ECSLs. 
The process is divided into three categories as defined below. 

A- The effect of the feedback strength variations: The effect of the increasing 
of the 77 in the laser dynamics could be divided into the two sections of low regime 77 
(such as rj < 0.0605) and high regime 77 (such as i] > 0.0727) [15]. In the low i] as it was 
expected | E | 2 dynamics reveals chaotic (CH), periodic (P1,P2,...), and quasi periodic 
(QP) behaviors. By the same token in the high regime | E | 2 dynamics reveals a regular 
plus package (RPP), QP and PI behaviors (see Fig. 1). The basic role of the 77 in laser 
dynamics is to displace of the periodic behavior windows, which can be presented by 
the Lyapunov exponent spectrum and the bifurcation diagrams. 

B- The effect of the feedback phase variations: The optical feedback phase 
is of a particular cyclic nature. Starting from a certain initial status, a variation of the 
C p by 27r must turn back to its initial status [T5] . To illustrate the effect of C p and its 
sensitive influence on the output intensity, | E | 2 , the study should focus on the two, 
Low and high regimes. In the low regime of 77, the output intensity | E | 2 presents a PI 
and P2, the chaotic behavior and QP (see Fig. 2). The appearance of the regular pulse 
packages (RPP) and disappearing of both the chaotic behavior and higher period (P2, 
P3,...) is the natural function of the output intensity | E | 2 in the high regime of the C p 
(see Fig. 3). The previous studies predicted bifurcation up the mode 10 (such as PI, 
P2, ...). However, based on the result of the present study, bifurcation can not exceed 
7 mode (see Fig. 4). 

C- The effect of pumping current variations: The bifurcation diagram in 
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variations of pump current follows the threshold value with respect to the r/ and C p . 
Where by considering the low rj, as can be seen in Fig. 5a, the | E | 2 is initially stable 
with a PI which undergoes early cascades of period-doubling to chaos. Chaotic window 
in bifurcation diagram varies based on the C p , as shown in Fig. 5(a-c). 
In high 77, as pictured in Fig. 6a, the | E | 2 is initially ( at first ) in PI dynamics, which 
by increasing the P, it undergoes RPP windows. As shown in Fig. 6(a-c), threshold 
of RPP is directly depended to C p with a higher and minor growth rate. At low 77, 
variations of periodic window and the increase of the chaotic window are the results of 
the increase of the pump current depicted in Fig. 7(a-c). The obtained results confirmed 
by plotting the time and phase portrait too (Fig. 8 and Fig. 9). Fig. 8(a-b) and Fig. 
9(a-b), verify the increase of the chaotic window. Fig. 8(c-d) and Fig. 9(c-d), on the 
other hand confirm variations of periodic window. Also based on the result the increase 
of the RPP windows and the disappearance of the QP in high 77 is the result of the 
increasing of the pump current (see Fig. 10 and Fig. 11). The overall picture of | E | 2 
dynamics has been presented in Fig. 12. 

5. Summery and Outlook 

Previous studies demonstrated have not provided sufficient insists towards the nature 
of ECSLs dynamics. Thus, there is a need to theoretically clarify its dynamics [18J. 
As noted before, bifurcation diagram can not reveal the hidden aspects of Laser 
dynamics [301 02] ■ Therefore, employing the Lyapunov exponent spectrum method can 
cover this methodological deficiency. The obtained results shed some lights on control 
process of the ECSLs dynamics. As a conclusions, it can be stated that r\ shows its 
influence on Laser dynamics by selecting a method of transfer form stable status to 
unstable status [15] such as (Periodic, QP and chaotic) behavior or (RPP, QP and PI). 
Therefor 77 can be regarded as the most important factor in stability of the Laser. Where 
Cp acts as a selection method in the scenario. Pump current is the other important factor 
in Laser dynamics, an essential factor which can influence the ECSLs dynamics from 
practical point of view. The use of pump current for controlling the ECSLs dynamics 
can also be considered a new field for further investigations. 
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Figure Captions 

Fig.l. Development of the dynamics for | E | 2 (P = 0.8 and C p = 1) (a)- Bifurcation 
diagram (b)-Lyapunov exponent spectrum. 

Fig. 2. Illustration of the effect of the C p variation on | E | 2 (P = 0.8 and rj = 0.036) 
(a)- Bifurcation diagram (b)-Lyapunov exponent spectrum. 

Fig. 3. Illustration of the effect of the C p variation on | E | 2 (P = 0.8 and n = 0.090) 

(a) - Bifurcation diagram (b)-Lyapunov exponent spectrum. 

Fig. 4. A comparative description of various periodic modes(P = 0.8) (a)- r\ = 0.042, 

(b) - v = 0.0455, (c)- 7/ = 0.048. 

Fig.5. The analysis of | E | 2 dynamics^ = 0.0455), (a)- C p = -1, (b)- C p = 0, (c)- 
C p = 1. 

Fig.6. The analysis of | E | 2 dynamics(?? = 0.115), (a)- C p = —1, (b)- C p = 0, (c)- 
C p = l. 

Fig. 7. The role of pump currents as initial condition in periodic and chaotic windows 
(rj = 0.0455) (a),(b) P = 0.6; (c), (d) P = 0.7; (e), (f)- P = 1.4. 

Fig.8. Time series of | E \ 2 {r] = 0.0455)(a)- C p = -1 and P = 0.6 (QP), (b)- 
C p = -1 and P = 0.7 (Chaotic), (c)- C p = -2.835 and P = 0.7 (P4), (d)- C p = -1 and 
P = 1.4 (P2). 

Fig.9. Phase diagram of | E | 2 (77 = 0.0455) (a)- C p = -1 and P = 0.6 (QP), (b)- 
C p = -1 and P = 0.7 (Chaotic), (c)- C p = -2.835 and P = 0.7 (P4), (d)- C p = -1 and 
P = 1.4 (P2). 

Fig 10. The effect of P in increasing the RPP windows(?? = 0.042) (a)-P = 0.8 (b)- 
P = 1.4. 

Fig 11. Increasing RPP windows on | E | 2 (C p = 0, 77 = 0.09 and P = 0.8) (a)- 
Time series (b)- Phase diagram. 

Fig 12. The overall picture of | E | 2 dynamics. 
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